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LUMPED  GEOPOTENTIAL  HARMONICS  OF  ORDER  14, 

FROM  THE  ORBIT  OF  1967-llG 

by 

D.  G.  King-Hele 

SUMMARY 

The  satellite  1967-llG,  which  had  an  orbital  inclination  of  40°, 
passed  through  14th  order  resonance  with  the  Earth's  gravitational  field 
in  1974.  The  changes  in  its  orbital  inclination  at  resonance  have  been 
analysed  to  obtain  values  for  four  lumped  14th-order  harmonics  in  the 
geopotential,  with  accuracies  equivalent  to  about  5  cm  in  geoid  height. 
Analysis  of  the  eccentricity  was  also  attempted,  but  did  not  yield  useful 
results . 

As  no  previous  satellite  analysed  at  14th-order  resonance  has  had 
an  inclination  near  40° ,  the  results  have  proved  to  be  valuable  in  deter¬ 
mining  individual  14th-order  harmonics  in  the  geopotential. 
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1  INTRODUCTION 

A  satellite  is  resonant  with  respect  to  the  Earth's  gravitational  field  when  its 
track  over  the  Earth  repeats  after  a  number  of  revolutions,  perhaps  14  or  15.  The  analy¬ 
sis  of  such  orbits  is  now  well  established  as  an  accurate  method  for  determining  geo¬ 
potential  harmonics  of  a  particular  order,  such  as  15th  order ^ .  The  geopotential  is 
expressed  as  a  double  infinite  series  of  tesseral  harmonics  of  order  m  and  degree 
S.(>m),  and  an  orbit  experiencing  mth-order  resonance  is  sensitive  to  harmonics  of  order 
m  ,  and  to  a  lesser  extent  those  of  order  2m  .  Analysing  a  resonant  satellite  yields 
values  of  'lumped'  harmonics,  which  are  linear  sums  of  individual  harmonic  coefficients 
of  order  m  .  If  a  number  of  these  mth-order  resonant  orbits  at  different  inclinations 
are  analysed,  a  number  of  values  of  lumped  harmonics  are  obtained,  and  the  resulting 
equations  can  be  solved  for  individual  harmonic  coefficients. 

Harmonics  of  order  14  are  derived  by  analysing  orbits  which  experience  14th-order 

resonance,  ie  the  satellite's  ground  track  over  the  Earth  repeats  every  14  revolutions. 

The  orbits  most  suitable  for  analysis  are  those  which  contract  very  slowly  under  the 

influence  of  air  drag  and  therefore  remain  close  to  14th-order  resonance  for  several 

years.  Such  orbits  arise  by  chance  from  time  to  time,  but  the  variety  of  inclination  so 

far  available  has  been  inadequate.  In  a  previous  determination  of  individual  14th-order 

2 

harmonics  from  resonance  analysis  ,  good  results  were  obtained  only  for  inclinations  of 
50°,  67°,  70°,  81°  and  87°.  However,  orbits  at  such  high  inclinations  are  affected 
primarily  by  harmonics  of  degree  less  than  20,  and  lower-inclination  orbits  are  essential 
if  the  evaluation  of  high-degree  harmonics  is  to  be  improved. 

In  1975  we  asked  for  US  Navy  orbital  elements  for  the  satellite  1967-1 1G,  inclina¬ 
tion  40°,  which  passed  through  14th-order  resonance  in  1974.  (It  was  a  fragment  from 

the  French  satellite  Diademe  1.)  These  orbital  elements  were  not  analysed  at  the  time  of 

2 

the  previous  determination  of  14th-order  harmonics,  because  1967-1 1G  passed  through 
resonance  quite  rapidly  and  the  orbital  elements  were  at  rather  wide  time  intervals. 
Disappointingly,  no  subsequent  satellite  with  an  orbital  inclination  near  40°  has  passed 
more  slowly  through  14th-order  resonance:  so  the  orbital  inclination  and  eccentricity  of 
1967-1 1G  are  now  analysed  to  obtain  values  of  lumped  harmonics  for  use  in  a  new  evalua¬ 
tion  of  individual  14th-order  harmonics.  As  expected,  the  analysis  was  quite  difficult 
and  not  fully  satisfactory;  however,  the  results  have  proved  to  be  good  enough  to  give 

3 

I  better-defined  values  for  individual  14th-order  harmonics  of  high  degree  , 

2  THEORY 

In  this  section  the  standard  theory  is  repeated,  so  as  to  define  the  lumped 
coefficients  and  .  The  longitude-dependent  part  of  the  geopotential  can  be 

4  ™ 

written  in  normalised  form  at  an  exterior  point  (r,  0,  A)  as 
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where  r  is  the  distance  from  the  Earth's  centre,  0  is  co-latitude,  X  is  longitude 

o  2 

(positive  to  the  east),  u  is  the  gravitational  constant  for  the  Earth  (398600  km'^/s  ), 

i>m 

1. 

and  S , 


R  is  the  Earth's  equatorial  radius  (6378.1  km),  P™  (cos  6)  is  the  associated  Legendre 


function  of  order  m  and  degree  5,  ,  and  C„  and  S_  are  the  normalized  tesseral 

“  J.m  Urn 

harmonic  coefficients,  of  which  those  of  order  m  =  14  and  28  are  relevant  here.  The 
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normalizing  factor  N.  is  given  by 
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The  rate  of  change  of  inclination  i  caused  by  a  relevant  pair  of  coefficients, 

2 

and  S  ,  near  B:a  resonance  may  be  written 


£ 

IT  =  iirr  (f)  ^ilmp‘^epq(‘^  i  "  q<^)l]  . 


(3) 


where 


is  Allan's  normalized  inclination  function  ,  G 


is  a  function  of 
6 


imp  "  **  *  ’  “ipq 

eccentricity  e  for  which  explicit  forms  have  been  derived  by  Gooding'^,  denotes 

'real  part  of  and  j  =  .  The  resonance  angle  i  is  defined  by  the  equation 


a((<)  +  M)  +  6(S  -  v) 


(4) 


where  oj  is  the  argument  of  perigee,  M  the  mean  anomaly,  the  right  ascension  of 
the  node  and  u  the  sidereal  angle.  The  indices  y,  q,  k  and  p  in  equation  (3)  are 
integers,  with  y  taking  the  values  1,  2,  3  ....  and  q  the  values  0,  ±1,  ±2,  ....; 
the  equations  linking  i,  m,  k  and  p  are:  m=YB;  k=Ya“q»  2p=i-k.  In 
equation  (3)  there  is  a  factor  (1  -  e'^)^  on  the  right-hand  side,  which  is  taken  as  1. 

Here  6  =  14  and  a  =  I  ,  and  the  m-suffix  of  a  relevant  (C„  ,  S.  )  pair  is 

decided  by  the  choice  of  y  .  The  values  of  I  to  be  taken  must  be  such  that  !.  >  m 

and  (i  -  k)  is  even.  The  successive  coefficients  which  arise  (for  given  y  and  q)  may 
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be  grouped  into  a  lumped  harmonic,  written  as 


59 'k 
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where  Z  increases  in  steps  of  two  from  its  minimum  permissible  value  JIq  ,  and  the 
are  constant  coefficients,  with  =  1  ,  The  values  of  the  can  be  obtained  from 

equation  (3),  and  R.H.  Gooding  has  written  a  computer  program  PROF  for  their  evaluation. 

The  rate  of  change  of  eccentricity  e  caused  by  the  (£,  m)  harmonic  near  6:a 
resonance  can  be  written^ 
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Again  there  is  a  factor  (1  -  e'^)  on  the  right-hand  side,  which  is  taken  as  1. 

The  relative  importance  of  the  many  possible  terms  arising  in  equations  (3)  and  (6) 
is  indicated  in  section  3. 

3  PRECONCEPTIONS  OF  THE  TERMS  IMPORTANT  FOR  1967-116 
3. 1  Inclination 

The  main  terms  contributing  to  di/dt  in  equation  (3)  are  listed  in  Table  1,  with 
the  factor  n(R/a)''  (14  cosec  i  -  cot  i)  taken  out.  For  1967-1 1G  at  resonance  this 
factor  has  the  value  301  radians/day. 


Table  1 


The  main  terms  in  the  equation  for  di/dt 


°  ^5,)4.7{Sh’ 

1  ^ /a\  14  - 

2  *  \r/  14  -  cos  i  14,14,7 

|cj^^  sin  (♦  -  id)  -  cos  (♦  -  ii))| 

i  JLl  /  a  \  14  -  2  cos  1  - 

2  ®  \R/  14  -  cos  i  14,14,6 

sin  (♦  +  «)-  cos  (♦  ♦  <>i)| 


<l)  ^28,28,13  (^8'  cos  2*} 

1A»  /iV*  ^  5 

\aj  14  -  cos  i  '29,28,14 

js^g*  sin  (2*  -  (d)  +  Cjg'  cos  (29  -  ii))| 
26-3  cos  i  - 

Va/  14  -  cos  i  '29,28,13 

|s^g’^  sin  (2*  *  u)  *  cjg’^  cos  (2*  ♦  u)| 


239  2  14  ■*■  cos  i  - 

8  ®  14  -  cos  i  'l5,14,8 


sin  (♦  -  2<,j)  ♦  C^*"’  cos  (♦  -  2u)| 


,  133  2  14  -  3  cos  i  5 

8  ®  14  -  cos  1  'l5,14,6 


sin  (♦  +  2u)  ♦  C"^’^  cos  (♦  +  2u)| 


The  orders  of  magnitude  of  the  lumped  coefficients  (C,  S)^!  in  Table  1  can  be 

A  •  '  ^  ”5  2 

assessed  very  approximately  by  assuming  that  C.  is  of  order  10  It  and  taking 

X,  ,m 


to  be  of  order 


•  ”Q  k 

similarly  for  S^’  .  Sometimes,  of  course,  the  real  values  of  the  lumped  coefficients 


may  be  small  by  chance,  and  this  method  will  overestimate  their  important.  Table  2  gives 
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the  values  calculated  in  this  way,  all  x  10 


Table  2 


Orders  of  magnitude  of  terms  in  Table  1,  for  1967-1 1G 


N.  Y 

q 

1 

2 

0 

7.1 

0.25 

1 

2.4 

0.2 

-1 

2.7 

0.2 

2 

0.5 

- 

-2 

0.5 

- 

From  Table  2  it  appears  that  only  the  terms  in  (y,  q)  =  (1,  0),  (1,  1)  and  (1,  -1)  are 
likely  to  be  of  importance. 

3.2  Eccentricity 

The  main  terms  contributing  to  de/dt  in  equation  (6)  are  listed  in  Table  3,  with 
the  factor  n(R/a)'  taken  out.  For  1967-1 1G  at  resonance  this  factor  has  the  value 
14.6  day  ' . 


2 


From  Table  A  it  appears  that,  for  e  ,  the  terms  in  (y,  q)  =  (1,  O  and  (1,  -1) 
will  be  the  most  important,  though  terms  in  (y,  q)  =  (1,  2)  and  (1,  -2)  may  also  be 
significant . 

An  alternative  method  for  evaluating  the  orders  of  magnitude  of  the  terms  would  be 

to  take  the  larger  of  the  C.  or  S.  in  a  particular  Earth  model.  Unfortunately, 

Jem  Jem 

with  1967-1 1G,  the  largest  are  for  2A  <  H  <  61  ,  where  the  values  in  existing  Earth 

models  are  not  considered  very  reliable. 

4  THE  PASSAGE  THROUGH  RESONANCE 

At  the  time  of  its  passage  through  resonance,  1967-1 1G  had  a  perigee  height  near 
540  km,  an  apogee  height  near  1100  km  and  an  orbital  period  of  101.3  minutes. 

There  were  59  sets  of  US  Navy  orbital  elements  available,  at  dates  between  19  April 
1974  (modified  Julian  day  42156)  and  10  January  1976  (MJD  42787).  During  this  time  the 
value  of  4>  ,  given  by  equation  (4),  changed  from  -3.99  deg/day  to  +8.68  deg/day,  with 
exact  resonance  ($  =  0)  on  9  September  1974  (MJD  42299).  Fig  1  gives  the  variation  of 
■1'  as  a  broken  line,  and  it  can  be  seen  that  its  rate  of  change  became  much  slower  in 
1975,  because  the  air  density  was  then  lower  than  in  1974,  as  a  result  of  lower  solar 
activity.  The  unbroken  line  in  Fig  1  shows  the  variation  of  $  over  the  same  time 
interval . 

Fig  1  shows  that  the  time  interval  covered  by  the  orbital  data  is  far  from 
symmetrical  about  the  exact  resonance,  and  this  has  several  implications  for  the  analysis. 
First,  the  terms  in  (y,  q)  =  (1,  -1),  that  is  the  terms  in  (I*  +  u) ,  are  likely  to 

have  an  appreciable  effect  only  near  the  time  when  ($  +  u)  is  nearly  constant,  that  is 

•  •  • 

when  $  =:  -  u)  .  But  w  =  6.3  deg/day  for  1967-1 1G,  so  the  main  effects  of  the  (y,  q)  = 

(1,  -1)  terms  occur  when  =  -  6.3  deg/day  -  well  before  the  start  of  the  orbital  data. 
Thus  it  is  likely  that  the  coefficients  of  the  (y,  q)  =  (1,  -1)  terms  will  prove  to  be 
indeterminate;  and  this  is  just  what  happens.  The  same  applies,  a  fortiovi,  with  the 
(y.  q)  =  (li  “2)  terms. 

So  the  terms  to  be  considered  in  the  fittings  are  (y,  q)  =  (1,0), (1,  1)  and  possibly 
(2,  0)  for  inclination,  and  (y,  q)  =  (1,1)  and  (1,2)  for  eccentricity. 

5  ANALYSIS  OF  INCLINATION 
5 . 1  Treatment  of  the  data 

The  US  Navy  orbital  elements  were  first  converted,  with  the  ELTRAN2  computer 

program,  to  midnight  epochs  in  the  five-card  format  of  the  RAE  PROP  orbital  elements. 

The  liinisolar  and  zonal  harmonic  pertubations  to  i  were  calculated  with  the  aid  of  the 

program  PROD^,  with  a  one-day  integration  interval  and  restarts  at  intervals  of  about 

40  days.  The  perturbations  indicated  by  PROD,  which  never  exceeded  0.0042“,  were  then 

subtracted  from  the  raw  values.  The  maximum  di/dt  due  to  solar  radiation  pressure 

perturbations  was  estimated  as  3  '<  10  ^  deg/day,  which  is  negligible  since  di/dt  due 

-4 

to  resonance  is  2  *  10  deg/day.  Perturbations  due  to  tides  were  also  ignored. 
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Fig  2  shows  the  values  of  inclination  cleared  of  perturbations.  Initially  all 
the  values  of  i  were  assumed  to  have  errors  of  0.003°, 

5.2  The  fittings  and  their  results 

These  values  were  then  fitted  with  an  integrated  form  of  the  theoretical  equation 

g 

(3) ,  with  the  aid  of  the  THROE  computer  program  ,  for  various  pairs  of  values  of 

(y,  q) ,  and  assuming  an  atmospheric  rotation  rate  A  =  0.9  rev/day.  The  first  run  was 

with  (y,  q)  =  (1,  0),  since  Table  2  indicates  that  these  should  be  the  dominant  terms. 

The  results  were  quite  good,  with  the  measure  of  fit  e  having  the  value  1.5  (we  define 

2 

c"  as  the  sum  of  the  squares  of  the  weighted  residuals,  divided  by  the  number  of  degrees 
of  freedom).  However,  it  was  obvious  that  some  of  the  values  were  ill-fitting  and  would 
remain  so  -  particularly  the  high  values  near  MJD  42600  visible  in  Fig  2.  The  accuracy 
of  a  number  of  these  values  was  relaxed  by  a  factor  of  2,  and  e  then  decreased  to 
0.96,  Next  the  (Yj  q)  =  (1,  1)  terms  were  added  and  slightly  improved  the  fit  (e  =  0.95) 
without  much  altering  the  values  for  the  (y,  q)  =  (1,  0)  terms.  Further  relaxations  in 
accuracy  were  then  made,  by  doubling  the  sd  of  any  values  with  residual  exceeding  2e  : 
this  reduced  e  to  0.80.  Finally  the  (y,  q)  =  (2,  0)  terms  were  added:  this  reduced  e 
to  0.75  and,  although  the  values  for  the  (2,  0)  terms  were  unlikely  to  be  reliable,  the 
improvement  in  fitting  seemed  worthwhile,  so  the  (2,  0)  terms  were  retained.  Some 
further  relaxations  and  restorations  in  accuracy  were  now  needed,  and  in  the  final  run 
the  assumed  sd  of  each  value  is  as  shown  in  Fig  2.  The  value  of  e  was  0.68,  and  the 
values  of  the  lumped  harmonics  which  emerged  are  as  follows: 


10^C°>’  = 

lU 

-65  ± 

163  , 

9-0  1 

10  =  -55  ±  164  , 

(7) 

io^c!’°  = 

1 4 

1.3  ± 

1.6  , 

10^s!;°  =  2.5  ±  1.3  , 

14 

(8) 

lO^C^'^  = 

10  1.23 

1 .4  ± 

0.4  , 

3-0  2 

10  =  -1.1  ±  0.5  . 

(9) 

The  standard  deviations  in  equations  (7)  to  (9)  are  rather  meaningless  as  they 
stand,  because  of  the  large  Q  values  of  1967-1 1G.  To  assess  the  significance  of  the 
standard  deviations,  it  is  worth  trying  to  give  an  indication  of  the  error  in  geoid 
height  to  which  they  correspond;  and  this  can  perhaps  best  be  done  by  dividing  each  sd 
by  the  largest  of  the  relevant  Q  factors  and  then  multiplying  by  the  Earth  radius 
(6378  km).  The  approximate  geoid  height  accuracies  for  the  three  pairs  of  standard 
deviations  in  equations  (7)  to  (9)  are  then  found  to  be  7  cm,  5  cm,  and  70  cm  respec¬ 
tively.  The  28th-order  harmonics  have  too  large  an  error  to  be  of  interest,  and  will  be 
ignored;  but  the  first  two  pairs  have  accuracies  similar  to  those  achieved  with  other 
orbits  analysed  at  14th-order  resonance. 

On  the  basis  of  the  10  ^ rule,  the  values  of  the  (C,  S)^’^  would  be 

-9  .  ^ 

600  >'10  ,  so  the  values  (7)  are  unusually  small  -  and  that  is  why  they  are  smaller 

than  their  sd.  The  10  ^ rule  gives  (C,  S)|?*^  as  of  order  6  *  10  so  the  values 


10 


0  2 

(8)  are  of  about  the  expected  magnitude.  The  (C,  S)  ’  terms  should  be  of  order 
-3  -5  2 

0.1  x  10  on  the  10  /Z  rule;  probably,  therefore,  the  values  (9)  are  considerably 
too  large. 

5 . 3  Equations  for  the  individual  harmonic  coefficients 

Each  of  the  lumped  harmonics  (C,  and  (C,  can  be  expressed  as  a 

linear  sum  of  individual  harmonic  coefficients,  as  indicated  in  equation  (5).  For 
1967-1 IG  the  explicit  forms  of  these  equations  are  as  follows 


?0.1.  c 

^14  •  ^ 


15,14  -  '■■“Cl7,14  *  '“■”S9,14  -  '5-»«21,14  *  '3-'“23,14  '  ‘•-5'25,14 


'■^^^27,14  *  ''■^'^29,14  "  ^*^^^31,14  "  '•^^^33, 14  *  ^■^^‘'35, 14 


=  (-65  ±  163)  T*  10 


(10) 


the  equation  is  the  same  as  (10)  with  C  replaced  by  S  and  the  numerical 
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value  on  the  right-hand  side  changed  to  (-55  ±  164)  x  10 


(11) 


^14  • 


S4,14  -  ’0-^S6,14  *  ^3.9C,g^,,  -  106.8C2o^,,  ^  '  ’78-5^24, 14 


+  103. 7C,,  ,,  +  8.6C,q  ,,  -  72.9C-.  ,,  +  51.1C-,  ,,  +  10.5C-.  ,, 
26,14  28,14  30,14  32,14  34,14 


-  39.4Cgg  =  (1.3  ±  1.6)  X  10 
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(12) 


the  equation  is  the  same  as  (12)  with  C  replaced  by  S  and  the  numerical 
value  on  the  right-hand  side  changed  to  (2.5  +  1.3)  x  lO  ^  ,  (13) 


The  series  on  the  left-hand  sides  are  terminated  at  degree  36  because  it  is  unlikely 
that  higher-degree  terms  can  be  successfully  evaluated  until  we  have  much  more  data. 

Equations  (10)  to  (13)  show  that  the  individual  coefficients  of  degree  19-24  are 
likely  to  make  the  largest  contribution  to  the  lumped  harmonics. 

6  /XNALYSIS  OF  ECCENTRICITY 

The  prospects  of  successfully  analysing  the  variation  in  eccentricity  seemed  poor 
from  the  outset,  because  there  were  such  large  changes  in  e  due  (presumably)  to 
perturbations  other  than  resonance.  The  values  of  e  after  correction  for  the  effects 
of  zonal  harmonics  and  air  drag  still  exhibited  a  largo  variation  -  from  near  0.039 
initially  to  more  than  0.040,  This  corresponds  to  a  change  in  perigee  distance  of  more 
than  7  km,  which  is  much  greater  than  would  be  expected  as  a  result  of  resonance.  In 
the  initial  THROE  run,  the  first  15  values  of  e  were  ill-fitting,  and  there  seemed  no 

harm  in  omitting  them,  because  the  resonant  variation  of  eccentricity  is  'centred'  on 

«  • 

the  point  where  1  =  w  ,  which  is  quite  late  in  the  data  set,  at  MJD  42550  (see  Fig  1). 


’0058 


1 1 

In  subsequent  shortened  runs  the  early  values  still  fitted  badly:  eventually,  only  33 
values  of  eccentricity  were  retained,  and  the  fitting,  with  (y,  q)  =  (1,  1),  is  shown  in 
Fig  3.  This  gave 

=  4.0  ±  1.7  ,  =  19.3  ±  1.7  ,  (i4) 

14  14 

with  £  =  1.16  ,  after  accuracies  were  relaxed  by  a  factor  of  2  for  all  values  having 
residuals  greater  than  2t  .  The  density  scale  height  was  taken  as  75  km. 

The  value  of  ^14^  expected  from  the  10  rule  is  6  x  10  ^  so  the  high  value 

in  equation  (14)  is  suspect,  and  the  likely  culprit  is  solar  radiation  pressure.  The 

largest  contribution  to  the  equation  for  the  variation  in  e  produced  by  solar  radiation 
9 

pressure  is  from  the  term  in  sin(n  +  u  -  L) ,  where  L  is  the  solar  longitude.  For 
1967-110,  the  angle  (11  +  w  -  L)  increases  at  a  rate  of  only  0.32®  per  day,  from  120®  at 
MJD  42479  to  220°  at  MJD  42787;  consequently  the  variation  in  e  due  to  solar  radiation 
pressure  closely  mimics  in  form  the  curve  of  Fig  3,  being  sinusoidal  with  a  maximum  near 
MJD  42670.  It  is  not  possible  to  remove  the  srp  effects  reliably,  because 

a.  the  mass,  size  and  shape  of  1967-1 IG  are  unknown; 

b.  its  decay  rate  is  irregular,  suggesting  irregular  variations  in  the  effective 
cross-sectional  area; 

c.  the  srp  cross-section  is  probably  also  variable,  and  may  vary  differently  to 
the  drag  cross-section;  and 

d.  the  srp  effect  is  probably  larger  than  the  resonance  effect  and  mimics  it  in 
form. 

2 

The  average  decay  rate  of  1967-1 IG  indicates  a  mass/area  ratio  of  about  20  kg/m  .  Taking 
this  value  for  the  srp  mass/area  and  using  only  the  main  term  in  the  srp  perturbation 
equation,  we  find  the  values  of  e  after  removal  of  srp  effects  are  as  shown  in  Fig  4. 
The  THROE  fitting,  shown  by  the  curve  in  Fig  4,  now  gives 

io^cl;“  =  -0.7  ±  1.6  ,  io^s!;°  =  7.5  +  1.6  ,  (15) 

14  *  14 

with  t:  =  1.10  .  The  values  (15)  are  quite  different  to  (14),  and  could  be 

reduced  to  a  value  consistent  with  (8)  if  the  srp  effect  was  increased  by  25%,  which  is 
about  the  probable  limit  of  error.  We  must  conclude  that  the  values  (14)  are  quite 
spurious,  because  of  the  neglect  of  solar  radiation  pressure,  while  the  values  (15)  need 
to  be  given  larger  sds  -  perhaps  about  ±4  -  to  take  account  of  the  possible  errors  in 
the  removal  of  srp  effects. 

In  Fig  4  there  is  an  oscillation  in  unison  with  u  ,  so  a  fitting  with  (y,  q)  = 

=  (1,  I)  and  (0,  1)  was  also  made,  and  is  shown  in  Fig  5.  This  gave 

10^c’’°  =  0.3  ±  1.2  ,  10^sl;°  =  9.7  ±  1.0  ,  (16) 

14  *14 
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with  e  =  0.80  ,  after  several  of  the  relaxed  sds  were  restored  to  the  standard  value. 
Though  nominally  more  accurate  than  (15),  the  sds  in  (16)  still  need  to  be  augmented  to 
allow  for  the  possible  mis-estimation  of  the  srp  effects,  and  again  values  of  about  ±4 
may  be  appropriate. 

All  in  all,  the  values  of  equation  (8)  seem  likely  to  be  more  reliable  than  (15) 
or  (16). 

7  CONCLUSIONS 

The  satellite  1967-1 1G  passed  through  14th-order  resonance  with  the  Earth's  gravi¬ 
tational  field  in  September  1974.  Analysis  of  the  changes  in  inclination  between  April 
1974  and  January  1976,  as  given  by  59  US  Navy  orbits  (Fig  2),  has  yielded  values  for 
four  lumped  harmonics  of  order  14,  and  the  accuracy  of  the  evaluation  is  equivalent  to 
about  5  cm  in  geoid  height:  see  equations  (7)  and  (8).  The  eccentricity  was  also 
analysed,  but  the  results  are  regarded  as  unreliable  because  of  uncertainties  in  the 
evaluation  of  perturbations  due  to  solar  radiation  pressure,  which  are  larger  than  the 
resonance  effects. 

Among  satellites  previously  analysed  as  they  passed  through  14th-order  resonance, 
none  has  had  an  inclination  between  34®  and  48®,  so  the  results  from  1967-1 IG  at  40® 
inclination  provide  new  information  which  should  be  useful  in  determining  individual 
harmonic  coefficients  of  order  14  -  particularly  those  of  high  degree.  This  expectation 

3 

has  proved  valid  . 
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